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Consistent statistical physical description is given for systems where the elementary excitations 
are composite objects. Explicit calculational scheme is constructed for the energy density and the 
total number of thermodynamical degrees of freedom, based on the spectral function of the system. 
One demonstrates through characteristic examples that single quasiparticle contributions combine 
non-linearly in these quantities. Relation to the Gibbs paradox is also discussed. 



QCD thermodynamics shows paradoxical faces in the 
crossover region connecting the hadronic and the quark 
phases. Variants of the perturbative approach, which 
provide a reasonable interpretation of the numerical sim- 
ulations in restricted temperature ranges, work with per- 
turbed ideal gas systems where the mean free path (^) 
is very large. On the other hand, as for example the in- 
terpretation of recent RHIC measurements suggest [1, 2], 
hadronic matter near Tc behaves as an almost ideal fluid 
where ^ — 0. Enforcing a free gas description to an ideal 
fluid, however, is possible only with strong interactions. 

But what seems to be strongly interacting in one way, 
can be weakly interacting in terms of the adequate de- 
grees of freedom. In PT one uses resummation to find 
them; and indeed, perturbative series are improved con- 
siderably eg. for entropy with HTL resummation [3], 
thermodynamics with screened PT [4] or with 2PI re- 
summation [5]. In effective approaches like the hadron 
resonance gas model or massive gluon picture, also non- 
interacting particles can give an account for complicated 
behavior in QCD like charm hadronization [6], equation 
of state [7-9], or interaction measure [10]. The success 
of these approaches suggests that probably most part of 
the strong interactions between free particles is incorpo- 
rated in the spectrum, and the residual interactions of the 
so-dressed excitations are small. 

We therefore may assume that the adequate descrip- 
tion is based on a model which contains weakly interact- 
ing excitations with nontrivial spectrum (not necessarily 
quasiparticles) . In this paper we consider the most simple 
effective theory described by a single real bosonic excita- 
tion, and omit all interactions. We primarily interested 
in the questions of how finite lifetime or multiparticle 
(threshold) effects can modify the thermodynamics, and 
how can different degrees of freedom show up in or vanish 
from the statistical ensemble dynamically. 

The treatment of an effective field theory with a quite 
general spectrum, however, requires to rethink some rel- 
evant questions concerning the foundations of statistical 
physics. The problem is the following: the partition func- 
tion 
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the multiplicities gn (for simplicity we assume no con- 
served charge, ie. no chemical potential here). In a 
fundamental theory the basic excitations are elementary, 
therefore their multiplicities are fixed permanently, while 
in an effective theory the thermodynamic degrees of free- 
dom might change dynamically. To have a cleaner char- 
acterization of the multiplicities, and also to approach 
the question of appearance/disappearance of modes in 
the ensemble, let us define the total number of (inter- 
nal) degrees of freedom as the sum of degeneracies in a 
1-particle case at zero spatial momentum 
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The value of Ndof is always an integer in the fundamen- 
tal theory and it appears directly in physical quantities: 
for example the pressure of bosonic gases at high tem- 
perature is P = NdofPsB, where Psb is the Stefan- 
Boltzmann (SB) limit Psb ^ TrT'^/90. 

In an effective theory, however, elementary excitations 
are in fact multiparticle states from the point of view of 
the underlying fundamental theory. Fig. 1/a shows three 
characteristic cases how the multiparticle nature mani- 
fests itself in the spectral function g. In case A the two 
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requires the knowledge of the energy levels En as well as 
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FIG. 1. a.) spectral function for: A.) well separated peaks, 
B.) merging peaks and C.) threshold. Figure b.) shows the 
behavior of the degrees of freedom from (11) as we go from 
case A to case B (and beyond). 

well-separated peaks each represent independent quasi- 
particles: they can be treated as a fundamental particles, 
and so Ndof ~ 2 here. In case B we should avoid the 
double-counting of the states in the overlap region, so we 
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expect Ndof < 2. If the widths F oo, finahy we get one 
peak and we expect N^of — >■ 1- The behavior of N^of as a 
function of the widths should be something like Fig. 1 /b 
- which results actually from an explicit evaluation (cf. 
below (11)). Case C is even more difficult, since there 
is no quasiparticle-like excitation at all; but we have to 
determine also its contribution to the thermodynamics. 

The problem of determining the correct degeneracies 
is in close connection with the old problem of indistin- 
guishability of particles in gases and the Gibbs paradox. 
Consider two energy levels with separation AE, then 
Niof = 2 for AE > 0, even for AE 0, but N^of = 1 
if AE = (since then wc have only one level). There- 
fore Ndof, and so the SB limit of the energy density is a 
non-analytic function of AE. However, according to the 
argumentation of the previous paragraph, for real gases 
with finite (eg. thermal) width, for AE < T the num- 
ber of degrees of freedom Ndof goes smoothly to one. A 
similar conclusion was drawn in [11]. By a better descrip- 
tion of the degeneracies in effective theories we may hope 
to have a deeper understanding of the Gibbs paradox in 
general. 

The key observation for the correct description is that 
statistical physics is determined by the dynamics, so we 
shall start with the Lagrangian of the effective theory. In 
case of a quadratic model the spectral function fixes com- 
pletely the dynamics, and so the thermodynamics and 
also Ndof. Below the spectral function will be the input 
of the analysis and the consequence of a few physically 
appealing non-trivial choices will be explored. 

The most general quadratic, spacetime translation in- 
variant action reads 



= \ j d'^xd'^y^p{x)K{x-y)ip{y). 
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The kernel K is to be constructed to represent the given 

spectral function. If it is continuous, K should be time- 
nonlocal. This immediately poses the question of causal- 
ity; but causality in this system simply means that the 
spectral function is zero for spacelike separation. This re- 
quirement is granted by construction, since the spectral 
function serves as input here. 

We should note that because of time-nonlocality some 
standard techniques will not work. First of all we do not 
have canonical formalism (or at least we would need in- 
finitely many fields). Also there is no direct connection 
between the imaginary time formalism and the thermo- 
dynamics: indeed, the standard derivation is based on 
the fact that the Hamiltonian density depends on the 
canonical momentum as 11^/2 [12]. In order to define 
thermodynamics we have to define the energy through 
the Noether current belonging to the time translation 
invariance, and then we should compute its expectation 
value at finite temperatures. We note that the formal 
correspondence between the time translation e~'^* and 
statistical operator e~^^ is still valid, and so KMS con- 



dition remains true: this will help us to take expectation 
values. Once we know the temperature dependent energy 
E{T), we can define the complete thermodynamics from 
that. 

To connect the kernel to the spectral function we com- 
pute the analytic continuation of the retarded/advanced 
propagator from the spectral function by the Kramers- 
Kronig relation 
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The kernel in Fourier space (denoted by /C), can be ob- 
tained from the relation: Gb(p) = IC~^{pQ + ie,p). Using 
the hermiticity of the action (3), there remains 



/C(p)= ReG-^\p). 
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If the spectral function is relativistic invariant, then the 
same holds for Gr and for the kernel, too. Prom now 
on we assiimc relativistic invariance (although for the 
applicability of the approach it is not necessary) . 

The next step is the construction of the energy- 
momentum tensor: it is the conserved current belonging 
to the spacetime translation symmetry. Using standard 
techniques [13] we obtain 
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where 

D^iu]C{p) = P/i^ - QnulC = 2pi^p^^ - g^^K. (7) 

The latter, relativistic form is explicitly symmetric. 

Expectation value of this form can be taken using the 
KMS condition iG^^{p) = n{po)e{p), where n is the Bose- 
Einstein distribution. We find 



(T^,) = J D^^Kip) Q +n(Po)) e{p), 



(8) 



where = ^ J . This is already position in- 
dependent. 

The leading 1/2 represents the vacuum energy. To 
obtain a finite expression we subtract the expectation 
value at zero temperature. 
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where 5g{p) = g{p) — Qo{p), and is the zero tempera- 
ture spectral function. 

In a rotationally invariant system only the diagonal el- 
ements survive. The 00 component is the energy density 



+ 



e = I DK{p) 
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n{po)g{p) + ^Sg{p) 
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where DJC{p) = DoolC{p). 

We note that s is independent on the normalization of 
the spectral function (if g Zq then Gr ZGr and 
K. K,/Z, so Z drops out). This means that only the 
energy levels count, not the normalization (which is the 
consequence of the definition of the fields). 

Having e(T) we can obtain the free energy from the 
relation e = the pressure from p = —f and the 

entropy density from s = ^{e+p). Prom s we can define 
also a quantity which can be interpreted as the number 
of internal degrees of freedom. At high temperatures 
iT'iPo) ^ T/po, and the average energy per mode is T. 
Therefore the number of modes for a given momentum is 
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p^DICip)g{p). (11) 
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This is similar to the "equivalent photon number" intro- 
duced by Weizsacker and Williams. N^of is dimension- 
less, temperature independent quantity which, as we will 
see below, in case of the discrete spectrum yields indeed 
the number of energy levels. 

After we have derived the relevant formulae we can 
apply them for some particular spectral functions. First 
of all one easily verifies that for a relativistic free particle 
(with spectral function q{p) = 2tt sgnpoS{p'^ — m^)) one 
obtains one degree of freedom (iV^o/ = 1) from (11), and 
the energy density from (10) yields the standard formula 



ei(m,T) = ^ j dpon{po) pl^fi^ 
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Consider now two stable particles with masses mi, 2 
and wave function renormalizations Zi^2- Since the result 
is independent on the global normalization of g, we can 
choose Zi + Z2 = 1 and consider 

g{p) = (27r) sgnpo [^i'5(/ - mj) + Z^Sip^ - m^] . 

(13) 



This yields Gb.{p) = Zi{p^ 
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with Landau prescription, and IC{p) = {p^ — m\){p'^ — 
m\)/ {jp^ - rfi^) where = Z2m\ + Zim%. Then we find 

dK ^ p^- 2p^m^ + Z2m\ + Zim\ 



dp- 



{p^ — w?)"^ 



If wc evaluate the internal degrees of freedom, we obtain 
A^do/ = 2 independently on the details, in particular on 
the weights of the Dirac-deltas. If, however, mi = 7712 at 
the beginning, then (like in the previous case), we have 
Ado/ = 1- Therefore we indeed count the number of 
the particle species, which is discontinuous: this is the 
manifestation of the Gibbs paradox in our expression. 
Similarly, the energy density reads 

^ r£l('71l,T) +£i(m2,T) ifTOi7^TO2 .^g. 

l£i(m,T), if mi = 777-2 = TO, 



where ei comes from (12). This formula is also non- 
analytic for TO2 — >■ TOi, and independent on Zi,2- 

To have a deeper insight into the calculations, we re- 
mark that in the energy density (11) we have the com- 
bination DtCg. Since /C = on the mass shells, we have 
to consider the dK/dp^ term. Near the mass shell it has 
the following property: 
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It is this non-interchangeability of the limits which lies 
behind the Gibbs paradox. A physical interpretation of 
these limits is that a finite width gives finite resolution 
of the spectra; if the peaks are closer than their widths, 
than then one cannot resolve them, and we are in the first 
limiting case. If they are farther, than we see two peaks, 
and the second limiting procedure has to be applied. 

The above observation generalizes to any number [N) 
of Dirac-deltas. Let us take in general 

e(po)= V^i(27r)5(po-Wi), Gfl = V 

■~ ^ Po - LOi + le 

(17) 

The kernel is its inverse, the derivative of the kernel reads 
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Since at pq — >■ Wj for any i the retarded Greens function 
diverges, its inverse is zero: ie. lC{p) = on any mass 
shell. On the other hand, if .^j ^ 
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and so the Zi factors drop out from the result. The 

number of degrees of freedom from (11) is Ndof = N, 
and the energy density is = X]„ei(TO„,T), both are 
independent on the normalization. 

After the infinite lifetime examples we consider the 
case of realistic gases. First let us consider a single degree 
of freedom with finite lifetime, represented by a spectral 
function of Breit-Wigner type 
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where ujp = 



The corresponding retarded Greens 



function and the kernel are G^{p) = {po + iT)-^ — 
and }C{p) = p^ — mP — F^. We can easily compute the 
number of degrees of freedom from (11) to obtain N^of = 
1! This means that, independently on the lifetime of the 
excitations, we always have a single degree of freedom. 

For the energy density we evaluate (10), the result can 
be seen on Fig. 2/a. One can see that for F/m as large 



4 




0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 






2SB 






r=o 


1 -delta 






r=0.2 


50*lhr 





1 1.5 2 2.5 3 
T 



b.) 



FIG. 2. a.) Energy density coming from a single Lorentzian 
for r = (ie. infinite lifetime) and F/m = 0.2. SB limit 
denotes the Stefan-Boltzmann limit, b.) Energy density for 
two Lorentzians with mi — 1 and m2 = 2 masses and Fi = 
F2 = F width. 2SB means double of the SB limit. 1-delta is 
a single Dirac delta case. 50*thr shows the 50-times energy 
density coming from a threshold. 



as 0.2, the energy density is close to the one for a stable 
particle. This observation makes possible to neglect the 
width for independent real quasiparticles in thcrniody- 
namical calculations. 

Although a single Lorentzian yields Ndof — 1, this does 
not mean - because of nonlinearity of the expressions (10) 
and (11) - that for N Lorentzians we obtain N degrees of 
freedom. To see this, let us take the linear combination 
of two such Lorentzians: 



Q2{p) = Zigi{p,Ti) + Z2gi{p,T2). 



(21) 



Then, for Ndof we obtain Fig. 1/b which we also dis- 
cussed in the introduction. If the width is small as com- 
pared to the distance of the peaks then we get Ndof = 2: 
in fact the F ^ limit is the case of two Dirac-deltas 
seen above. If the width gets larger, the two Lorentzians 
start to overlap, and we continuously arrive to the case 
of one broad Lorentzian seen above with Ndof = 1- The 
phenomena seen on Fig. 1/b. is nothing else than the 
Gibbs-paradox, smeared out for real gases. 

The energy density can be seen on Fig. 2/b. What 
we see is that for mi — 1 and TO2 = 2 a particle width 
r = 0.2 is already very close to the 1-particle case. We 
therefore see the reduction of degrees of freedom also in 
the energy density. 

Finally we can take a pure threshold with g{p) ~ 
Y^l — rn? I]p- . With m 1 we see the result for the 
energy density on Fig. 2/b. In order to see something 
at all, we multiplied the result by 50 - this suggests that 
thresholds have very tiny contribution to the energy den- 



sity. This result can be interpreted as the fate of multiple 
bound states: when they are melted to a broad contin- 
uum, then their contribution vanishes from the EoS. 

As a conclusion we stress that one has to be cautious in 
treating effective quasiparticle models, even in the weakly 
interacting case. Since quasiparticles are composite ob- 
jects, they do not represent independent degrees of free- 
dom, only in the limit when their mass difference is much 
larger than their width. As a consequence their contri- 
bution to different physical quantities like energy density 
do not simply add up: we observe nonlinear behavior, 
when two finite width quasiparticles can melt into a sin- 
gle one, or multiparticle states can give vanishing con- 
tribution (cf. Fig. 2). We expect that this nonlinear 
behavior may play important role also in other physical 
effects, like charm suppression or transport phenomena 
in fluids. Our method can also be used to identify other 
currents and compute transport coefficients from the ef- 
fective theory. This is the subject of ongoing projects. 
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